One Health

Background: Previous studies indicated that regional fluctuations in COVID-19 impact were associated with socioeconomic, demographic and environmental factors early in the pandemic. It is unclear whether factors associated with regional COVID-19 mortality and SARS-CoV-2 infections changed over the course of the pandemic. Methods: We evaluated the association of ecological factors with municipality-specific COVID-19 mortality in the national population of Austria between January 2020 and June 2023. Analyses were conducted separately per calendar year and pooled across years via random-effects meta-analysis. Outcomes were analysed using univariable and multivariable zero-inflated negative binomial (ZINB) regression, including annually averaged environmental exposome exposure (including air pollution and UV radiation), socioeconomic, demographic and infection dissemination predictors with significance threshold set to p<0.005. Secondary analysis modelled SARS-CoV-2 infections using of negative-binomial model (NB). Two additional approaches for smaller predictor set selection were applied as sensitivity analyses. Results: Within the entire observation period 18,430 COVID-19 deaths were reported. In univariable analyses (pooled RR (99.5% CI)), COVID-19 mortality was significantly associated with the proportion of commuters (0.79 (0.73-0.86)), residents aged 75-84 (1.18 (1.01-1.37)) and [&ge;]85 years (1.29 (1.04-1.60)), residents with at most compulsory education (1.15 (1.02-1.29)), employed residents aged 15-64 (0.89 (0.85-0.93)), and homeowners (0.91 (0.83-0.99)), as well as with precipitation (0.94 (0.89-0.99)) and the average time with a mean temperature of 5-15 C (0.95 (0.90-0.99)). In multivariable analyses the share of population aged 85 or older (1.30 (1.09-1.54)) and home ownership (0.92 (0.87-0.98)) were independently significant. Mortality predictor estimates were relatively stable across pandemic years. Across 5,963,882 documented SARS-CoV-2 infections share of employed in 15- to 64-year-olds (1.01 (1.01-1.02)), part-time employment rate (1.01 (1.00-1.02)) and district level nursing care allowance recipiency (0.97 (0.95-0.98)) were associated with detected infections. Overnight stays, population 85 or older, income, compulsory education or less and air pollution all had multiple significant estimates with opposite directionality, leading to non-significant overall associations. Sensitivity analyses outcomes mirrored main analyses outcomes. Conclusions: Throughout the pandemic in Austria, demographic and socioeconomic factors were associated with COVID-19 mortality and infections, but many associations, especially for infections, were unstable over time. These findings underscore the importance of structural vulnerabilities in epidemic impact.

Background: Lean Healthcare has been adopted as an approach to reducing waste, enhancing process efficiency, and optimizing clinical and administrative workflows. However, empirical evaluations aimed at establishing a causal relationship between Lean Healthcare interventions and outcomes remain limited. This evidence gap is particularly pronounced in the context of large-scale implementations in low- and middle-income countries, where more efficient allocation could yield substantial improvements in population health. Methods: We conducted an impact evaluation study employing a staggered difference-in-differences approach to estimate the effects of the Lean Healthcare in Emergency Departments (ED) initiative in the Brazilian context. The intervention was implemented between January 2017 and December 2025 across 309 hospitals. We analyzed hospital-level indicators of the number of hospitalizations and the average length of stay (in days) among patients admitted through urgent or emergency pathways. Results: Our estimates suggest that the intervention was associated with an average increase of 38.6 admissions per hospital per month (p = 0.033), corresponding to approximately a 5.9% increase relative to the monthly admission volume. We also observed a reduction of 0.184 days in the average length of stay (p = 0.058), equivalent to an approximately 3.0% decrease in hospitalization time. These effects exhibited distinct temporal dynamics: the effect on the number of hospitalizations emerged gradually and increased with the duration of the intervention, whereas the reduction in length of stay became evident earlier. Conclusion: This study provides evidence of a causal effect of Lean Healthcare implementation in emergency departments on operational performance. The simultaneous increase in admissions and reduction in length of stay suggests improved bed productivity and effective capacity without expanding physical infrastructure. These findings support the use of flow-management methodologies as an effective strategy for mitigating emergency department overcrowding while increasing service capacity.

Importance: Coccidioidomycosis, also known as Valley fever, is a serious fungal disease caused by inhaling Coccidioides spores aerosolized from the soil. Despite rising public health concern, the distribution and incidence of coccidioidomycosis in the US remain poorly characterized as not all states are required to report cases, and case under-ascertainment is presumed to be high, particularly outside of the endemic region. Objective: To estimate monthly, county-level coccidioidomycosis incidence across the US while accounting for widespread under-ascertainment. Design, Setting, and Participants: We integrated multisource data, including human cases diagnosed in electronic health records, serological testing of dogs, and ecological and sociodemographic predictors, in a machine learning model estimating monthly, county-level coccidioidomycosis incidence. The model was trained on surveillance data from Arizona and California from 2013 to 2023, then applied to predict incidence in states across the contiguous US over the same time period. Main Outcomes and Measures: Incidence of coccidioidomycosis at the monthly county level. Results: The model predicted coccidioidomycosis incidence from 2013 to 2023 for 3,108 counties from all states in the contiguous US. The model identified 60 counties in non-reporting states where predicted incidence was [&ge;]2 cases per 100,000 population (considered moderate to high incidence). In states with mandatory reporting outside of Arizona and California, the model predicted county-level incidence rates that were 28% higher (95% bootstrap interval, 0.1%-69.8%) than those captured in national surveillance, equating to a total of 608 potentially missed cases per year (95% bootstrap interval, 340-964). Conclusions and Relevance: In this cross-sectional study, the model quantified case under-ascertainment for coccidioidomycosis relative to Arizona and California, and highlighted the potential for alternative data sources to fill in gaps in mandatory reporting. Estimating incidence while accounting for under-ascertainment and non-reporting may help efforts to improve public and physician awareness and to detect emergence in new regions.

Background. Nigeria's reported Lassa fever case fatality ratio (CFR) rose from 18.9% at week 23 of 2025 to 25.0% in 2026, prompting concerns over failing clinical care. Because surveillance artefacts can also alter CFRs, candidate mechanisms that published data can resolve were identified and evaluated. Methods. This ecological study extracted weekly state-level counts from 474 Nigeria Centre for Disease Control and Prevention situation reports (2017-2026). Per a pre-specified plan, years were compared at phase-matched points, year-on-year CFR changes were decomposed into within-state and between-state components, and tipping-point ascertainment shifts calculated. Weekly incident and cumulative published totals were compared post hoc. Findings. Of 11 candidate mechanisms, 2 could be evaluated fully and 5 in part. On weekly incident counts, 2026 did not differ from 2025 at the phase-matched midpoint (+0.3 percentage points, 95% CI -5.0 to 7.6); within-state (-1.1, 95% CI -4.7 to 3.0) and between-state (+0.5, -2.0 to 2.9) components spanned zero. Conversely, the published cumulative CFR was 5.4 to 6.2 percentage points higher because 70 deaths were added retrospectively by week 36 outside weekly reports. Reconciling cumulative figures would require 279 missed survivors (6% of suspected cases not confirmed). Interpretation. The 2026 CFR rise came from deaths added to cumulative totals outside the weekly counts, and published data cannot attribute it to care. Endemic zoonotic CFRs cannot serve as clinical outcome measures unless testing volumes, deaths added late, and post-mortem sampling are published alongside case counts. Funding. US NSF-NIH-NIFA EEID; UK BBSRC; Birgitta Sintring Fellowship (Uppsala University).

This paper examines the impact of Cancer Care Ontario's regionalization of lung cancer surgery on blood transfusion use for thoracic surgery at The Ottawa Hospital. Using aggregate data scraped from the paper ''Changing trends in blood transfusion: an analysis of 244,013 hospitalizations'' by Shehata et al. (2014) we estimated the effect regionalization on the log mean number of blood transfusions among thoracic surgeons (who treat lung cancer), utilizing other surgical specialties as controls. We found that the average number of blood transfusions increased significantly among thoracic surgeons following the opening of a new cancer care centre at the hospital in December 2009, while control specialties showed no parallel change. Contextualizing these results in a modified ''population level'' treatment threshold model, it is argued that the main driver of this increase in blood transfusions among thoracic surgeons was due to the influx of less skilled thoracic surgeons to The Ottawa Hospital, and the closing of thoracic surgery units at smaller hospitals which was part of the broader regionalization process.

Introduction: Pesticide exposure has been associated with poorer mental health, yet little is known about whether seasonal fluctuations in agricultural pesticide use are associated with cyclical changes in mood among people living in agricultural communities. We evaluated associations between 10 peak and off-peak pesticide spray seasons and ecological momentary assessment (EMA)-derived mood in a floricultural region in Ecuador. Methods: 496 participants completed repeated EMA surveys across 10 assessment administrations over 2.5 years (2022-2024), including one single-assessment administration and nine 10-day bursts (Nobservations=14,332). Mood outcomes included negative affect (sadness, anxiety, stress, anger, tiredness, difficulty concentrating; overall mean=2.79, Standard Deviation [SD]=1.32), positive affect (happiness, motivation; overall mean=4.21, SD=1.60), and behavioral dysregulation (aggression and impulsivity, overall mean=1.86, SD=1.29). Generalized linear mixed models estimated pooled, between-person and within-person associations between peak pesticide spray seasons and mood, adjusting for demographic and longitudinal factors Effect modification by sex, floricultural occupation, and cohabitation with a floricultural worker was assessed. Results: Peak spray seasons were associated with higher negative affect ({beta}=0.11 [0.06, 0.16]) and behavioral dysregulation ({beta}=0.10 [0.05, 0.15]) but not positive affect ({beta}=0.01 [-0.05, 0.07]). Associations were driven by within-person changes. Sex modified associations with negative affect, which were stronger in males than females. Floricultural occupation modified associations, with poorer mood, greater behavioral dysregulation and lower positive affect among floricultural workers than non-workers. Discussion: Peak pesticide spray seasons were associated with worsening of negative affect and behavioral dysregulation. EMA may provide a sensitive approach for detecting short-term neurobehavioral changes associated with seasonal pesticide exposure.

Wastewater-based epidemiology uses both targeted and non-targeted bioanalytical approaches to detect and track disease-causing microbes in wastewater. Targeted approaches such as quantitative polymerase chain reaction (qPCR) remain widely used while sequencing-based approaches have slowly gained popularity. However, concordance between established targeted molecular assays and sequencing approaches remains incompletely characterized. To assess the validity of sequencing as an alternative to qPCR, we compared outcomes of both methods targeting five respiratory viruses in wastewater samples collected across 16 Texas cities over 4 years (May 2022-July 2026) of continuous sampling. We found positive correlations between qPCR and sequencing measurements, with both methods independently correlating with statewide clinical emergency department visits. qPCR positivity also increased with genomic sequencing coverage, and sequencing demonstrated greater potential than qPCR as an early respiratory virus outbreak detection method. Our findings support metagenomic sequencing as a sensitive and highly specific alternative to targeted qPCR for wastewater surveillance, and reinforce its utility towards long-term longitudinal monitoring of respiratory viruses.

Background: Both over- and under-antibiotic use pose major global health challenges, yet there is no consensus on how to quantify deviations from optimal prescribing or set evidence-based targets. We developed a benchmarking tool for inpatient antibiotic use and evaluated its application in 10 countries. Methods: We developed decision-tree models to estimate expected empiric antibiotic use for 11 indications, based on the WHO AWaRe book for two age groups: 28 days to 12 years and >12 years. We parameterised the models using country-specific antimicrobial resistance estimates in the 10 countries and cross-sectional Global Point Prevalence Survey (Global-PPS) case-mix data from 2021 for 10 countries and from 2015-2022 in a subset of five countries. We then compared estimated expected antibiotic use with observed empiric use (from Global-PPS) in defined daily doses for adults and days of therapy for children, across indications and AWaRe categories. Findings: In 2021, among patients aged >12 years, in all countries but one, observed empiric Access antibiotic use, as a proportion of total antibiotic use, was lower than expected and empiric Watch antibiotic use was higher than expected. Similar patterns were observed in all countries but one in patients aged 28 days to 12 years. The magnitude of discrepancies between expected and observed Access and Watch antibiotic use varied across indications and countries. The median observed-to-expected ratios of usage across countries ranged from 0.3-0.8 for Access antibiotics, and 0.7-5.3 for Watch antibiotics in patients aged >12 years in 2021. Temporal trends from repeated surveys from 2015 to 2022 showed little change in observed-to-expected ratios for Access and Watch antibiotics. Interpretation: This benchmarking tool uses Global-PPS data to quantify differences between observed antibiotic use and WHO AWaRe recommendations. Our tool can be adaptable to local guidelines, provides policy-relevant estimates of expected antibiotic use to support monitoring, procurement planning, and targeted stewardship intervention.

Vegetation plays a critical role in shaping fine-scale tick density by structuring local microclimates, habitat resources, and host interactions that largely determine tick survival and activity, with implications for tick-borne disease transmission. Tick burden in the United States continues to increase, with annual reported tick-borne disease cases more than doubling from 2004 through 2019. These trends call for a broader scope of tick-mitigation strategies that go beyond traditional prevention measures that focus on personal protection; and there is an opportunity to inform habitat-level interventions. Here, we quantify fine-scale vegetation structure across 48 subplots at tick surveillance sites in Western Massachusetts to assess local habitat features associated with Ixodes scapularis density and potential for human-tick contact risk. In assessing understory and overstory vegetation attributes, this study identified life stage-specific differences in predicted tick density, with higher adult densities in areas with deeper leaf litter and presence of invasive plant species and nymphs associated with high presence of ferns. Both adult and nymphal Ix. scapularis were predicted at higher densities in relatively low basal area forests (18.9 to 26.9 m2/hectare). This study design provides a simple, low-cost framework for organizations looking to incorporate tick-informed forest measurement ecology into their public health tick surveillance protocols with minimal training.

Tick-borne disease risk reflects the combined influence of vector abundance, pathogen prevalence, human exposure, and seasonal activity, yet these components are rarely integrated at the county level. We developed species-specific ecological risk scores for three medically important tick species (Acari: Ixodidae): Amblyomma americanum (L.), Dermacentor variabilis (Say), and Ixodes scapularis Say, across all 102 Illinois counties using data from 2019 to 2022. County-level predicted tick abundance and pathogen prevalence were scaled and combined using equal weights to generate primary vector-level risk scores. Robustness to the equal-weight assumption was evaluated by systematically varying the relative weights assigned to abundance and pathogen prevalence, while a separate sensitivity analysis incorporated reported tick-bite incidence where available. Associations with human disease were evaluated using negative binomial regression of tick-borne disease case counts on the corresponding risk scores, with county population included as an offset. Monthly phenology indicators were constructed by integrating human disease cases, reported tick bites, and sampled tick abundance and were evaluated in relation to temperature and precipitation. Seasonal risk was estimated by combining county-level vector risk with species-specific monthly phenology, generating monthly risk surfaces that were deployed through an interactive R Shiny application. Higher risk scores were significantly associated with greater disease case rates for each tick species. County rankings were highly stable under moderate changes in abundance and prevalence weighting. Seasonal activity differed among species, with A. americanum and D. variabilis showing concentrated warm-season risk and I. scapularis exhibiting a broader seasonal profile. These county-level seasonal risk maps identify where and when relative tick-borne disease risk is elevated in Illinois, supporting targeted prevention, surveillance, and clinical awareness.

Background: Commercial spaceflight now carries civilians on short missions, yet their stress-axis physiology is poorly described. Earlier shuttle data suggest that sympathetic responses predominate after short flights and glucocorticoid responses after long ones. The neutrophil-to-lymphocyte ratio (NLR) is a low-cost readout of both. Methods: Preregistered secondary analysis of open clinical laboratory data (NASA Open Science Data Repository OSD-569 and OSD-575) from the four Inspiration4 crew members, each sampled seven times: three draws before launch, one the day after return (R+1) and three during recovery (R+45 to R+194). Log-transformed NLR was modelled with a linear mixed model with a random intercept per crew member and Kenward-Roger inference. Lymphocyte and eosinophil counts and C-reactive protein were secondary outcomes with false discovery rate correction. Results: All 28 draws were analysed. NLR at R+1 did not differ from pre-flight (fold change 1.04; 95% CI 0.81 to 1.34; p = 0.73); two crew members increased (+25% and +46%) and two decreased (-27% and -11%). During recovery NLR was 23% lower than pre-flight (fold change 0.77; 95% CI 0.64 to 0.92; p = 0.005). No secondary outcome changed after correction; lymphocytes were 18% lower at R+1 (95% CI -35% to +4%; q = 0.31). Between-person differences explained 58% of the variance. An ordinary regression ignoring repeated measures missed the recovery effect (p = 0.08). Conclusions: A 3-day orbital flight produced no consistent leukocyte stress signature one day after return. The lower NLR during recovery suggests that pre-flight values, not recovery values, deviate from the resting state, which has implications for how baselines are defined in short-duration spaceflight studies.

Japanese encephalitis virus (JEV) is a leading cause of viral encephalitis in humans in Asia. Prior to 2021, only sporadic outbreaks had been detected in Australia's tropical northern extremes. However, during 2021-23, JEV infection was identified in 45 people and in more than 85 piggeries, together with many detections in feral pigs, mosquitoes, and other livestock across Australia, including in the temperate far south. Here we estimated the spatial extent of virus transmission during the outbreak (the realised niche) and the area environmentally suitable for transmission (the potential niche). We collated a national dataset of JEV surveillance and detection in humans, mosquitoes, domestic pigs, feral pigs, and sentinel chickens from 2021-23. We combined these data with (1) previously published predictions of habitat suitability for vertebrate host and vector species, and (2) estimated densities of humans, domestic pigs, and feral pigs, in a multipart assessment of spatial risk. Our model projected broad potential environmental suitability over much of northern and eastern Australia during the outbreak period, but floodplain wetland areas in the north and the inland east of Australia had both the highest estimated environmental suitability for transmission, and the strongest evidence for virus circulation.

Wolbachia-based interventions were less effective in achieving introgression in Rio de Janeiro than other settings. Reduced introgression occurred in locations where dengue was historically high, raising concerns of inequity. We examined factors affecting Wolbachia introgression in Rio through trajectory modeling, spatiotemporal Bayesian models, and Ae. aegypti abundance estimation. During the intervention (2017-2020), Wolbachia introgression was 26% (median; IQR,14-44%), with four distinct trajectories. Despite receiving more releases, densely populated areas with low socioeconomic status had lower introgression. In spatiotemporally explicit multivariable regressions, higher Ae. aegypti abundance, temperature, population density and lower socioeconomic status were associated with reduced introgression. Improved characterization of Ae. aegypti abundance and other structural/environmental predictors of introgression could be used to calibrate release doses to ensure equitable benefits from Wolbachia interventions.

Introduction: Acute gastroenteritis (AGE) is a cause of childhood morbidity and may be associated with human adenovirus (HAdV) infection. This study aimed to investigate the occurrence of HAdV in fecal samples from children with AGE in Salvador, Bahia, Brazil, as well as to molecularly characterize the identified strains and evaluate their cytopathic effects in cell culture. Methods: A total of 131 samples were analyzed using the nested polymerase chain reaction (Nested-PCR) technique. Molecular characterization was performed through Sanger sequencing and phylogenetic analysis, while Vero cell cultures were used to investigate the occurrence of cytopathic effects. Results: Of the samples analyzed, 13 (9.9%) tested positive for HAdV, while 118 (90.1%) were negative. The highest detection frequency was observed in children aged three years or younger, particularly those under one year of age. Phylogenetic analysis identified variants related to HAdV-F41 and HAdV-B3, revealing genetic differences and distinct clustering patterns compared with reference sequences. In cell culture, cytopathic effects were observed four days after inoculation, although further analyses are required to confirm productive viral replication. Conclusion: These findings demonstrate the presence of HAdV in the samples studied and contribute to understanding its genetic diversity and circulation in Salvador, highlighting the importance of epidemiological surveillance and further studies on its role in the etiology of childhood AGE.

Introduction: Raised blood pressure is an important cardiovascular risk factor in Cambodia; however, national averages may obscure substantial local geographic heterogeneity. We assessed the spatial distribution of raised blood pressure and its demographic, metabolic, behavioural, built-environment, and environmental correlates using the 2023 Cambodia STEPwise approach to noncommunicable disease (NCD) risk factor surveillance (STEPS) survey. Methods: We analysed 289 georeferenced primary sampling units (PSUs) from the 2023 Cambodia STEPS Survey. PSU-level case counts and numbers of adults with valid blood pressure data were modelled using a frequentist binomial generalised additive model with a two-dimensional Gaussian-process smooth. An explanatory model included selected PSU-level and contextual covariates; a separate mapping model used spatially continuous predictors to generate predictions on a 5-km grid. Design-based national and residence-specific prevalence estimates were calculated separately using the survey design and Step-2 weights. Results: Among 4,225 adults with valid blood pressure data, unweighted prevalence of raised blood pressure was 28.6%, compared with 18.8% after Step-2 weighting; weighted prevalence was 17.5% in rural areas and 21.7% in urban/suburban areas. Mean observed unweighted PSU-level prevalence was 29.0%. In the adjusted spatial explanatory model, the proportion aged 60-69 years ({beta}=1.969, SE=0.327; p<0.001) and overweight/obesity prevalence ({beta}=0.996, SE=0.265; p<0.001) were positive ecological correlates. The spatial model had lower Akaike information criterion and residual deviance than its non-spatial counterpart, although Bayesian information criterion was higher. Prediction uncertainty was greater in some sparsely sampled areas. Conclusions: Raised blood pressure showed significant geographic heterogeneity across Cambodia. Older age composition and overweight/obesity prevalence were positive PSU-level correlates, not individual-level causal effects. Predicted prevalence should be interpreted alongside uncertainty and should complement, rather than replace, design-based STEPS prevalence estimates for geographically informed NCD planning.

Introduction: Acute gastroenteritis (AGE) is a global public health problem and a major cause of mortality in low- and lower-middle-income countries. Human bocavirus (HBoV) is among the viral agents associated with this condition. Therefore, investigating the presence of HBoV in cases of AGE in Salvador, Bahia, as well as its prevalence and circulating strains, is important. Objective: To investigate the occurrence of HBoV infection during an AGE outbreak in Salvador. Methods: A total of 178 samples were analyzed for HBoV using nested polymerase chain reaction (nested PCR), and the amplification products were visualized by agarose gel electrophoresis. HBoV was detected in 5.7% of the samples analyzed. Conclusion: This study contributes to the understanding of HBoV epidemiology as a secondary agent associated with acute gastroenteritis during a norovirus outbreak in Salvador. Continued investigation and identification of the etiological agents involved in outbreaks are essential for evaluating the causes of AGE. These findings also highlight the need for broader diagnostic approaches and the importance of phylogenetic analysis for monitoring and controlling viral infections.

Wildfire smoke poses a growing health risk worldwide, yet the magnitude and temporal evolution of mortality after episodic exposure remain poorly resolved, limiting burden assessment. Using monthly county-level wildfire PM2.5 and mortality data from the United State, we exploit episodic exposure contrasts to deconvolve overlapping mortality responses and recover the impulse-response function. The excess mortality response peaks approximately three years after an exposure pulse and subsides by five years. We estimate 1.7 (95% uncertainty interval: 1.1-2.3) million smoke-associated deaths during 2011-2023, approximately ten times estimates using long-term ambient PM2.5 mortality relative risks. Circulatory and respiratory diseases contribute comparably and dominate this burden. Responses are markedly greater among adults [&ge;]65 years; during the COVID-19 pandemic, population-wide cumulative responses are significantly higher over at least the first seven months after exposure. Reducing wildfire smoke exposure can thus yield substantial health benefits over multiple years, particularly when population vulnerability is heightened.

Brazil has one of the highest dengue burdens globally, with transmission influenced by climate events. Yet, the effects of compound climate extremes on dengue risk and prediction have not been quantified to date. Using a Bayesian spatiotemporal model of monthly dengue cases across 435 health regions from 2010 to 2025, we estimated relative risk associated with interacting long-lag and short-lag drought indicators and temperature anomalies across nine Koppen-Geiger climatic zones, adjusting for absolute temperature and El Nino-Southern Oscillation. Long-lag dry, mid-lag warm and short-lag wet conditions were associated with the highest relative risk in five climatic zones, with temperature anomalies modifying hydrometeorological dynamics more strongly in temperate than tropical areas. Climate covariates improved the prediction of case counts and outbreaks compared to baseline models at lead times of 1, 3, 6 and 12 months, and for the upcoming season 5 months ahead. Improvements were greatest in the Northeast, Southeast, and Midwest. Recent dengue incidence provided the greatest improvement at the 1-month lead, whereas combining incidence with climate information performed best at 3- and 6-month leads. We extended our framework to weekly state-level predictions of the season ahead to align with national public health planning. We proposed a multi-staged dengue forecasting approach to enhance early warning and response across Brazil.

Background: While health system costs of COVID-19 are well documented, evidence on out-of-pocket costs to patients and families, including non-medical and time burden, is limited. This study aimed to estimate out-of-pocket medical, non-medical, and time burden of COVID-19 illness directly from patients and caregivers. Methods: Two cross-sectional surveys were fielded to national nonprobability samples of U.S. adults in Fall 2023: a personal burden survey (adult sample) and a caregiver burden survey (child sample). Surveys collected data on care setting, visit duration including travel and wait time, and medical and non-medical out-of-pocket costs. After applying quality-based exclusion criteria, 409 adult sample and 52 child sample responses were analyzed. Results: Over 87% of respondents were insured. Median out-of-pocket medical costs per visit were low, ranging from $5 (emergency department) to $10 (urgent care) in the adult sample, and from $0 (emergency department) to $75 (hospitalization) in the child sample. Median non-medical costs were $28 or less across care settings. Median time spent receiving care ranged from 90 minutes (primary care, child sample) to 300 minutes (emergency department, adult sample). Despite the relatively low reported out-of-pocket costs, COVID-related financial challenges were reported by 34% of the adult sample and 29% of the child sample, with the most common difficulty being paying non-medical bills (11% and 15%, respectively). Twelve percent in both groups reported problems paying COVID-19-related medical bills. Conclusions: In this post-pandemic non-probability sample, individual out-of-pocket costs for COVID-19 illness were relatively modest. However, time burden was substantial and one-third of respondents reported financial hardship. These estimates, including non-medical and time burden not captured in administrative data, can support more robust economic evaluations of COVID-19 prevention and treatment interventions.

Background: Cost-effectiveness analyses (CEA) built on simulation models are usually delivered through project-specific scripts. Each new model, and often each new question, rebuilds the same scaffolding: parameter management, calibration, base-case analysis, sensitivity analysis and export. This is costly, and is where traceability between an estimate and the model version producing it is lost. We present THALASSA, a web application providing this infrastructure once, for any model satisfying a small programming contract. Methods: A simulation model is connected by one file implementing a fixed contract of R functions that determine its behaviour. THALASSA exposes the complete CEA workflow in a browser: a generated overview of the connected model, parameter editing, calibration, base-case comparison with the efficiency curve, and deterministic and probabilistic sensitivity analysis at a willingness-to-pay threshold, each with timestamped file export carrying the parameter values it was run with. Every session is bound to a repository version, so results carry their generating model version. Results: We demonstrate THALASSA on two public models of contrasting difficulty: a simple tutorial model and a generic cancer screening model with two competing natural-history pathways and a 21-dimensional calibration deliberately built to be more complex. Both ran on the same unmodified setup, implementing the same contract, with no code changes. Conclusions: Separating analysis infrastructure from model logic makes the CEA workflow reusable across projects and reproducible by construction. A modelling team can move from a versioned repository to audit-ready, exportable CEA outputs without writing analysis code. The integration effort does not grow with model size and one installation serves unrelated projects.